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Proton Spin Content From Lattice QCD* 
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We calculate the form factor of the quark energy momentum tensor and thereby extract the quark orbital 
angular momentum of the nucleon. The calculation is done on a quenched 16^ x 24 lattice at /3 = 6.0 and with 
Wilson fermions at «: = 0.148, 0.152, 0.154 and 0.155. We calculate the disconnected insertion stochastically 
which employs the Z2 noise with an unbiased subtraction. This proves to be an efhcient method of reduce the 
error from the noise. We find that the total quark contribution to the proton spin is 0.29 ± 0.07. From this we 
deduce that the quark orbital angular momentum is 0.17 ± 0.08 and predict the gluon spin to be 0.21 ± 0.07, i.e. 
about 40% of the proton spin is due to the glue. 



(N 

in 
o 
00 
o 
a^ 



X 



To understand the spin content of the proton 
remains a challenging problem in QCD Ex- 
perimental [|| , 1^ and lattice results Q , ^ sug- 
gest that the quark contribution (^S) to the pro- 
ton spin is about 25 ± 10%. But to date, we have 
very little knowledge about the remaining part 
of the proton spin. We do not have reliable esti- 
mate about the spin contribution from the gluons 
or the orbital angular momentum of the quark. In 
this talk, I will show our lattice results on the to- 
tal angular momentum of the quarks and thereby 
deduce the quark orbital angular momentum and 
predict the gluon contribution to the proton spin. 

Recently it was shown Q| that one can decom- 
pose the total angular momentum of QCD in a 
gauge invariant way, i.e. 



d^xip^x X {-iD)}ip 



+ / d^x[xx {E X B)]. (1) 
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The forward matrix element of this operator in 
the proton defines the decomposition of the pro- 
ton spin i = -is + + J, 
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where is the 



quark spin contribution, Lq is the quark orbital 
angular momentum and Jg is the total angular 
momentum of the glue. To calculate the angu- 
lar momentum of the quark and the gluon in the 



proton, one first notices that the gauge invariant 
quark-gluon energy momentum tensor B is 



rpflU 



q ~^ 9 



i [ij^^''iD'''>ij + tp-f^'^iD'^^i;] 



(2) 



where the first part is the quark energy momen- 
tum tensor and the second one is that of the 
gluon. Form factors of this energy momentum 
tensor current can be defined as 

+ r2'?^s(g2)p(/^ia'')"q„/2m 

+ -Ts'''{q^Wq-^g''-q') 
m 

+ T,'^^^{q^)g^^m]u{p',s% (3) 



where p^' = {pf" + p'^)/2, = 



p'^ and u{p) 



*Talk presented by N. Mathur at Lattice '99, Pisa, Italy 



is the nucleon spinor. It is proved |g] that the 
total angular momentum of the quark or gluon is 
related to the sum of the Ti and T2 form factors 
at zero momentum transfer, i. e. 

J ^ <P,s\\ e-^^ J d^x (Tf^x^ - TO/gX^) \p, s > 
'■^ < P, s\p, s > 

= ^ [J'f (0) + r|'^(o)] = r9'ff(o). (4) 
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Therefore, by calculating the form factor T'^'3{q^) 
at different and then extrapolating to ^ 
limit one can obtain the total angular momentum 
of the quark and the gluon separately. In this talk 
I shall present our calculation of the total angular 
momentum of the quark, Jq = + Lg. 

To compute the matrix element, one needs 
to calculate the two and three point functions. 
Three point function has two parts: connected 
insertion (CI) due to valence and cloud quarks, 
and disconnected insertion (DI) arising out of sea 
quarks. For CI, we evaluate the form factor T{q^) 
in terms of the following ratio (for notation see 
ref. §): 



Tr 


^mGNTOj N{t2, tl, 0, —q) 


Tr 


^eGNN{tl,0) 


Tr 




Tr [reGNNiti,q)' 



^^e.kmqkT^^'iq'). (5) 

For CI, we calculate T'^{q^) at different q^ and 
then extrapolate that to the — * limit. Results 
are obtained for relatively light quarks with k = 
0.148, 0.152, 0.154 and 0.155. We follow the same 
technique used in [Q and Q. Fig.l shows the 
dipole fitting of T(q^ at different q^ for k = 0.152 
and 0.155. Plots for the other k's are similar. 
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Figure 1. Dipole fitting for CI. Extrapolated 
{q^ — » 0) values are shown by Tl{{)). 



The chiral limit for J^^ is taken with a lin- 
ear dependence on the quark mass m^a for these 
four K, (fig. 3). To account for the correlations, 
both fittings in q^ and niq are done by using the 
covariant matrix and the final error at the chiral 
limit is obtained from the Jackknife procedure. 
Finally, to get the continuum result wc multiply 



the lattice Tl(0) by the tadpole improved renor- 
malization constant for the operator T^-' and 
obtain the CI part of the total angular momen- 
tum of quarks J^^ = 0.43 ± 0.07, which almost 
saturates the proton spin. The calculation is done 
on a quenched 16^ x 24 lattice at /3 = 6.0 with 
Wilson fermions and with 100 configurations. For 
DI contribution, we use the relation 

1 V- rf°G'p^^^p(t/T) tf»a r 

^1^2^ — F^^^^^TT^ — ^ c + . Tj,, (6) 

i = l T Gpp{tf) 

where C is a constant and DI is obtained from 
the slope of the sum on the left hand side as a 
function of i/ . To calculate the trace of the three 
point function in DI, we follow the same method- 
ology as in ref. Q and use the same stochastic 
algorithm, which employs the Z2 noise estimator 
and provides the estimate with minimum vari- 
ance. In addition, we use two more techniques 
to reduce the errors. First one is the removal of 
the unwanted real or imaginary part of the three 
point function through charge conjugation and 
hcrmiticity (CH) symmetry. Secondly, we use an 
unbiased subtraction method ||l^ to reduce the 
variance due to the stochastic estimator, i. e. we 
use the following estimator 

p 

Tr {A-')=E[< ryt {A'' 0« )v >] , (7) 

where ?7's are Z2 noise vectors, O^'^'s are a set 
of P traceless matrices and Ai's are variational 
constants which need to be tuned to reduce the 
variance of Tr(yl.^^). In principle, one can take 
any set of traceless matrices as a choice of O^'^K 
But we choose a set of traceless matrices obtained 
from the hoping parameter expansion of the prop- 
agator in order to match the off-diagonal behav- 
ior of the matrix so that they can cancel the 
off-diagonal contributions to the variance [|o|. A 
variational program which minimizes the variance 
from the 100 gauge configurations is used to de- 
termine the optimum set of A's. With these two 
techniques, we obtain a reduction of error by as 
much as 3-4 times. The subtraction method is 
also proved to be very efficient for the calculations 
of fcrmion determinants lflO| and for the traces of 
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many other operators Fig. 2 shows a plot of 
sum vs tf (Eq.@) for k — 0.154 and = 1 with 
and without subtraction. Finally, we fit Tj^^ as a 
function of by a monopole form which is shown 
in Fig. 2 for k — 0.154. Similar to CI, we also 
use covariant matrix fitting and the final error- 
bar is obtained by the jackknife method. A finite 
mass correction factor is introduced while ex- 
trapolating to the chiral limit. For strange quark 
contribution we follow the same procedure as in 
ref. a. 
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Figure 2. (a) The ratio of Eq.(|§) for k = 0.154 
at — 1. Ratios with subtraction are shifted 
a little towards right, m is the slope, (b) The 
monopole fitting of T{q^) 
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Figure 3. The lattice T^j{0) and r^j(O) for CI 
and DI as a function of quark mass ma. The 
chiral result is indicated by •. 



From DI calculation, we find that J^^ , like 
the DI part of the quark spin h^,^^ , is also 
flavor symmetric, i.e., = J^^ ~ J^^ = 

-0.047 ±0.013. Therefore, the total DI contribu- 
tion J^^ = -0.14 ± 0.04. Here one should note 
that in ref. [g, results for iE^^ = -0.18 ± 0.03, 
which suggests that the sea quarks give very little 
orbital angular momentum contribution. Adding 
CI and DI contributions, we obtain Jg — 0.29 ± 



0. 07 and therefore we can predict the gluon an- 
gular momentum content from the spin sum rule, 

1. e., Jg = l/2-0.29±0.07= 0.21±0.07. From our 
results one can draw the following conclusions: 

1. The total angular momentum content of the 
quark in the proton is Jq = 0.29 ± 0.07, i.e., 
about 60% of the proton spin is contributed by 
the quarks. Since the lattice result in ref. Q gives 
E = 0.25 ± 0.12, one can deduce that the quark 
orbital angular momentum is 0.17 ± 0.08. Thus, 
one concludes that about 25% of the proton spin 
comes from the quark spin and about 35% comes 
from the quark orbital angular momentum. 

2. Gluon angular momentum contribution is pre- 
dicted to be Jg = 0.21 ± 0.07, i.e., about ~ 40% 
of the proton spin is attributable to the glue. 

3. For the sea quarks, almost all the contribution 
comes from the quark spin and the total sea quark 
contribution nearly cancels the contribution due 
to glue. As a result, the proton spin is almost 
saturated by the CI (valence & cloud quark con- 
tribution) of the quark angular momentum alone. 
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